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Abstract. Astrophysical black hole candidates are thought to be the Kerr black hole pre- 
dicted by General Relativity. However, in order to confirm the Kerr-nature of these objects, 
we need to probe the geometry of the space-time around them and check that observations 
are consistent with the predictions of the Kerr metric. That can be achieved, for instance, 
by studying the properties of the electromagnetic radiation emitted by the gas in the ac- 
cretion disk. The high-frequency quasi-periodic oscillations observed in the X-ray flux of 
some stellar-mass black hole candidates might do the job. As the frequencies of these oscil- 
lations depend only very weakly on the observed X-ray flux, it is thought they are mainly 
determined by the metric of the space-time. In this paper, I consider the resonance models 
proposed by Abramowicz and Kluzniak and I extend previous results to the case of non-Kerr 
space-times. The emerging picture is more complicated than the one around a Kerr black 
hole and there is a larger number of possible combinations between different modes. I then 
compare the bounds inferred from the twin peak high-frequency quasi-periodic oscillations 
observed in three micro-quasars (GRO J1655-40, XTE J1550-564, and GRS 1915+105) with 
the measurements from the continuum-fitting method of the same objects. For Kerr black 
holes, the two approaches do not provide consistent results. In a non-Kerr geometry, this 
conflict may be solved if the observed quasi-periodic oscillations are produced by the reso- 
nance vq : v r = 3 : 1, where vq and v r are the two epicyclic frequencies. It is at least worth 
mentioning that the deformation from the Kerr solution required by observations would be 
consistent with the one suggested in another recent work discussing the possibility that steady 
jets are powered by the spin of these compact objects. 
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1 Introduction 

In 4-dimensional General Relativity, uncharged black holes (BHs) are described by the Kerr 
solution, which is completely specified by two quantities: the mass, M, and the spin param- 
eter, a* = J/M 2 , where J is the spin angular momentum, of the compact object 1 [1—3]. A 
fundamental limit for a Kerr BH is the bound |a*| < 1, which is the condition for the exis- 
tence of the event horizon. For |a*| > 1, there is no horizon and the Kerr metric describes 
the gravitational field around a naked singularity, which is forbidden by the Weak Cosmic 
Censorship Conjecture [4]. 

It is widely believed that the final product of the gravitational collapse is a Kerr BH, as 
follows from a body of analytical and numerical studies [4-8] (but see [9-12] and references 
therein). Astronomers have also discovered several good astrophysical candidates [13]. In 
particular, we have about twenty stellar-mass BH candidates in X-ray binary systems. For 
the time being, we can only get a robust measurement of the mass of these objects. That 
is possible by studying the orbital motion of the stellar companion. The latter is relatively 
far from the compact object and its trajectory is well described by Newtonian mechanics. 
In this way, we can estimate the mass of the compact object without any assumption about 
its nature. As these masses turn out to exceed 3 Mq (the maximum mass for a neutron 
or quark star assuming no phase transition to exotic forms of matter at densities below the 
nuclear one [14, 15]), these objects are thought to be Kerr BHs, as they cannot be explained 
otherwise without introducing new physics. 

In order to confirm the Kerr-nature of the known astrophysical BH candidates, we need 
to probe the geometry of the space-time around them and check if it is consistent with the 
predictions of General Relativity [16]. That can be achieved in a number of way. In the case 
of stellar-mass BH candidates, we may get information on the geometry of the space-time 
around these objects by studying the thermal spectrum of their accretion disk during the 
high-soft state [17, 18], by the analysis of relativistic lines [19, 20], and possibly even by 

1 Throughout the paper, I use units in which Gn = c = 1. 
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measuring the power of their jets [21, 22]. If the stellar companion is a radio pulsar, its 
orbital motion can be accurately tracked and we can get information about the nature of 
the BH candidate [23]. Ground-based gravitational wave detectors are supposed to be able 
to observe BH quasi-normal modes (QNMs) in a near future and, since the ones of a Kerr 
BH depend only on M and J, the possible detection of at least three modes can test the 
Kerr-nature of the source [24]. Other approaches can instead be used only to probe the 
geometry of the space-time around super-massive BH candidates in galactic nuclei [25-34] . 

Another promising tool to test the nature of stellar-mass BH candidates is represented 
by the so-called quasi-periodic oscillations (QPOs). The X-ray power density spectra of some 
low-mass X-ray binaries show some peaks; that is, there are QPOs in the X-ray flux 2 [36]. 
These features can be observed in systems with both BH candidates and neutron stars and 
the frequencies of these oscillations are in the range 0.1-10 3 Hz. High-frequency QPOs in 
BH candidates (~40-450 Hz) are particularly interesting, as they depend only very weakly 
on the observed X-ray flux and for this reason it is thought they are determined by the 
metric of the space-time rather than by the properties of the accretion flow. If that is 
correct, high-frequency QPOs can be used to probe the geometry around stellar-mass BH 
candidates. For the time being, however, the exact physical mechanism responsible for the 
production of the high-frequency QPOs is not known and several different scenarios have been 
proposed, including hot-spot models [37, 38], diskoseismology models [39-42], and resonance 
models [43-46]. In these models, the frequencies of the QPOs are directly related to the 
characteristic orbital frequencies of test-particles. 

The possibility of using QPOs to test the strong gravitational field around stellar-mass 
BH candidates has been already studied [47-49]. In [47, 48], the authors consider a very 
specific case, the braneworld BH metric proposed in [50], in which the geometry of the space- 
time is fairly similar to the one of an electrically charged BH of General Relativity. Let us 
notice, however, that this solution may be strongly constrained by current Solar System ex- 
periments (but only if the Birkhoff's Theorem holds). In Ref. [49], the authors study QPOs 
in a quasi-Kerr metric, valid only for small spin parameters and small deviations from the 
Kerr solution. In this paper, I study the phenomenon in a quite generic stationary and ax- 
isymmetric non-Kerr background and I find a much more complicated picture. In particular, 
the possibility of the existence of vertically unstable orbits, absent in the Kerr background, 
allows for a larger number of possible combinations between different modes. I then consider 
the scenario of resonance models proposed by Abramowicz and Kluzniak [43-46], mainly 
motivated by the observed 3:2 double peak high-frequency QPOs in micro-quasars, and I 
compare theoretical predictions with observational data. Lastly, I compare these results with 
the measurements inferred from the continuum-fitting method of the same objects. In the 
Kerr case, this approach is equivalent to the comparison of the spin measurements from 
high-frequency QPOs and continuum-fitting method and the question of consistency of the 
two techniques has been discussed in Ref. [51]. While high-frequency QPOs and continuum- 
fitting method do not seem to provide consistent results for Kerr BHs, the conflict may be 
solved if the space-time around BH candidates deviates from the Kerr solution. 

The content of the paper is as follows. In Sec. 2, I discuss the characteristic orbital 
frequencies of test-particles in the Johannsen-Psaltis (JP) metric [52], which is a stationary 
and axisymmetric space-time with very generic deviations from the Kerr background. In 
Sec. 3, I apply the results of the previous section to the resonance models. In Sec. 4, I consider 

2 QPOs have been observed even in the spectra of super-massive BH candidates, see e.g. Ref. [35]. 
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the observational data of the three stellar-mass BH candidates with observed twin peak high- 
frequency QPOs and with a measurement of the mass. These objects are GRO J1655-40, 
XTE J1550-564, and GRS 1915+105. Fortunately, for them we have also measurements of 
the thermal spectrum of their disk during the high-soft state and I can thus compare the 
allowed region in the plane spin parameter-deformation parameter inferred from the QPOs 
with the one determined from the continuum-fitting method. Discussion and conclusions are 
in Sec. 5. 

2 Characteristic orbital frequencies 

Let us now focus on equatorial orbits, as the accretion disk around BH candidates is normally 
expected to lie on the equatorial plane of the system. Circular orbits of test-particles are 
characterized by three frequencies: the Keplerian frequency z^k (which is the inverse of the 
orbital period), the radial epicyclic frequency u T (the frequency of radial oscillations around 
the mean orbit), and the vertical epicyclic frequency ug (the frequency of vertical oscillations 
around the mean orbit). These three frequencies depend on the geometry of the space-time 
and on the radius of the orbit. While they are defined as the characteristic frequencies of the 
orbital motion for a free particle, there is a direct relation between these frequencies and the 
ones of the oscillation modes of the fluid accretion flow. 

In Newtonian gravity with potential V = —M/r, the three characteristic frequencies 
has the same value: 

M . . 

v K = v r = v e = . (2.1) 

They are plotted as a function of the radial coordinate in the top left panel of Fig. 1. 

In General Relativity, a key-ingredient is the existence of the innermost stable circular 
orbit (ISCO). In the Schwarzschild metric, the ISCO radius is nsco = 6M. In the Kerr 
background, the ISCO radius decreases (increases) as the spin parameter a* increases for 
corotating (counterrotating) orbits. For a maximally rotating Kerr BH (a* = 1), the ISCO 
radius is respectively nsco = M for the corotating case and nsco = 9M for the counterro- 
tating one. Circular orbits with radii smaller than nsco are radially unstable. That means 
the radial epicyclic frequency v r reaches a maximum at some radius r max > nsco and then 
vanishes at the ISCO. Keplerian and vertical epicyclic frequencies are instead defined up to 
the innermost circular orbit (also called photon orbit). Circular orbits with radius smaller 
than the one of the photon orbit do not exist. In the Kerr metric, vq > u r , while, for corotat- 
ing orbits, uk > vq- &V> an d vq in the Schwarzschild space-time are shown in the top right 
panel of Fig. 1. The Kerr cases with spin parameter a* = 0.9 and 0.998 and for corotating 
orbits are shown respectively in the bottom left and bottom right panels of Fig. 1. 

In the case of a generic non-Kerr background, the picture is more complicated. The most 
important difference is that the ISCO radius may be determined by the stability of the orbit 
along the vertical (instead of the radial) direction [17, 53]. In some metrics, there is also the 
possibility of the existence of one or more regions of stable circular orbits inside the ISCO and 
therefore with smaller radii. These regions of stable orbits are separated by a gap from the 
"traditional" ISCO; the details depend on the specific metric, but usually they do not play 
any physical role, as they can accept only particles with energy higher than the one at the 
ISCO. In what follows, I will consider the JP metric [52], as it seems to include all the relevant 
features of a generic space-time deviating from the Kerr solution. While several authors have 
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Figure 1. Profiles of the Keplerian frequency vk, the radial epicyclic frequency v r , and the vertical 
epicyclic frequency vg in Newtonian gravity with potential —M/r (top left panel), in the Schwarzschild 
space-time (top right panel), and in the Kerr space-time with spin parameter a* = 0.9 (bottom left 
panel) and a* = 0.998 (bottom right panel). See text for details. 



proposed specific alternatives to Kerr BHs for the astrophysical BH candidates (like Q-stars, 
gravastars, or non-Kerr BHs solutions of specific and theoretically motivated extensions of 
General Relativity), here the choice of the JP metric with a single deformation parameter 
can be motivated as follows. The aim of this work is not to test a specific theoretical model 
and to determine/constrain its parameters, but to investigate in a model-independent way 
possible deviations from the Kerr geometry of the space-time around astrophysical BHs. The 
JP metric is used to perform a null-experiment: it is like the Kerr metric with a deformation 
parameter measuring possible deviations from the Kerr solution and the spirit is to determine 
this deformation parameter and check it is zero. Current and near-future data are indeed not 
so good to map the space-time around a BH candidate and a single deformation parameter 
is used to figure out if the gravitational force is stronger or weaker than the one around a 
Kerr BH with the same mass and spin. Actually, the typical situation is even worse and 
with one measurement we can only infer one parameter: if we assume the Kerr background, 
we find the spin a* (but we find a wrong value if the object is not a Kerr BH), if we have 
also a deformation parameter, we constrain some combination of a* and of the deformation 
parameter. With two independent measurements, we can determine both the spin and the 
deformation parameter. The choice of the JP metric should thus be understood with this 
spirit and it makes sense if it is used to check if the geometry of the space-time is described 
by the Kerr solution, not to investigate the actual nature of the compact object or the exact 
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deviations from the Kerr background. 

In Boyer-Lindquist coordinates, the JP metric is given by the line element 



ds 2 
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The JP metric has an infinite number of deformation parameters and the Kerr solution 
is recovered when all the deformation parameters are set to zero. However, in order to 
reproduce the correct Newtonian limit, we have to impose eo = ei = 0, while €2 is strongly 
constrained by Solar System experiments [52]. In this paper, I consider a simple case with 
a sole deformation parameter 63 and e, = for % 7^ 3. A different choice of the deformation 
parameter, like £4 or £5 instead of €3, would not change the qualitative features of our non- 
Kerr metric, as well as our results and conclusions [54]. 

In a generic stationary, axisymmetric, and asymptotically flat space-time, the charac- 
teristic orbital frequencies can be computed numerically as follows. The line element of the 
space-time can be written in the canonical form 



ds 2 = g tt dt 2 + g rr dr 2 + geedO 2 + 2g tl j ) dtd(j) + g^dcf 2 , 



(2.4) 



where the metric components are independent of the t and <f> coordinates, which implies the 
existence of two constants of motion: the conserved specific energy at infinity, E, and the 
conserved z-component of the specific angular momentum at infinity, L z . This fact allows to 
write the t- and ^-component of the 4- velocity of a test-particle as 



u 



Eg<l><f> + L z g t< p 



u • 



9tcj> ~ 9ttg<p<t> 

From the conservation of the rest-mass, g^u^u 1. 



Egg, + L z g tt 
9uf> ~ 9tt9cf>cf> 

-1, we can write 



g rr r 2 + g ee 2 = V eS (r, 9, E, L z ) , 



(2.5) 



(2.6) 



where r = u r = dr/dX, 
V e ff is given by 



it 



d9/dX, A is an affine parameter, and the effective potential 
E 2 g ( f }< f ) + 2EL z g t<j) + L 2 g tt 



9ttj> ~ 9tt9H 



1. 



(2.7) 



Circular orbits on the equatorial plane are located at the zeros and the turning points of the 
effective potential: r = 9 = 0, which implies V e s = 0, and r = 9 = 0, requiring respectively 
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d r V e g = and dgV c g = 0. From these conditions, one can obtain the Kepler ian angular 
velocity Q,k, E, and L z of the test-particle: 



dd> -d r g t( j> ± J {d T gt4>f ~ (d r g tt ) (d r „ 

° K = J = a? • (2i 



E 9tt + gt^ ^ 

\J-gtt - 1g t(j $l - g^ti 2 ' 
V ~5tt - - 



where in Qk the sign + is for corotating orbits and the sign — for counterrotating ones. The 
Keplerian frequency is simply v& = VLk/2it. The orbits are stable under small perturbations 
if d 2 V cS < and djV^ < 0. 

The radial and vertical epicyclic frequencies can be quickly computed by considering 
small perturbations around circular equatorial orbits, respectively along the radial and verti- 
cal direction. If 5 r and 8g are the small displacements around the mean orbit (i.e. r = vq + 8 r 
and 9 = tt/2 + Sg), we find they are governed by the following differential equations 

' /J ' V + n% = , (2.11) 



dt 2 
dt 2 



+ n 2 9 S e = 0, (2.12) 



where [45] 



tf- 1 f2 13) 

r ~ 2g rr {ut) 2 dr 2 ' {lA6) 
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The radial epicyclic frequency is thus u r = £l r /2-K and the vertical one is vg = Qq/2tt. The 
radial profile of the three characteristic orbital frequencies in the JP background with spin 
parameter a* = 0.8 and four different values of the deformation parameter €3 is shown in 
Fig. 2. 

3 Resonance models 

In four stellar-mass BH candidates, we observe two high-frequency QPOs. It is remarkable 
that in all the four cases the ratio of the two frequencies is 3:2, suggesting a strong correlation 
between them. Possible theoretical models should thus be able to explain this feature. The 
resonance models proposed by Abramowicz and Kluzniak [43-46] seem to be quite appealing 
from this point of view. In Eqs. (2.11) and (2.12), the radial and the vertical modes are 
decoupled. However, it is natural to expect that in a more realistic description there are non- 
linear effects coupling the two epicyclic modes. In this case, the equations can be written 
as 

d26r +n 2 5 r = n 2 F r (5 r ,5g, d 4, d 4) , (3.1) 



dt 2 ' r " "T-ryr,»a, & , ^ 

'•^ +n 2 5g = n 2 Fg(8 r ,5g, d ^,^) . Ci.2, 



dt 2 ' ^vyr^v, dt , dt 
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Figure 2. Profiles of the Keplerian frequency vk, the radial epicyclic frequency v r , and the vertical 
epicyclic frequency vg in the JP space-time with spin parameter a* = 0.8 and deformation parameter 
£3 = 2.0 (top left panel), e 3 = 1.0 (top right panel), £3 = 0.8 (bottom left panel), and e 3 = —2.0 
(bottom right panel). See text for details. 



where F r and Fg are some functions that depend on the specific properties of the accretion 
flow. If we knew the details of the physical mechanisms of the accretion process, we could 
write the explicit form of these two functions and solve the system. Unfortunately, that is 
not the case. The strategy is therefore to guess possible properties and consequences of these 
equations and see if they can be fitted in a plausible physical scenario. 

3.1 Parametric resonances 

A simple but interesting scenario is to imagine that vertical oscillations are governed by the 
Mathieu equation [44, 46]: 

^ + n 2 e s e = -njh cos(n r t)8 e , (3.3) 

which corresponds to the case with F r = and Fg = —5 r 5g: the solution of Eq. (3.1) is 
simply 5 r = hcos(£l r t) and so we obtain Eq. (3.3). The Mathieu equation describes indeed a 
parametric resonance with 

^ = -, n = 1,2,3,... (3.4) 
ilg n 

The resonance is stronger for smaller values of n. In the Kerr background, vg > u r , and 
therefore the resonance n = 3 can naturally explain the observed 3:2 ratio if the upper 
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frequency u\j is associated with vg and the lower frequency z/l with v r . In the JP space-time 
with €3 > 0, vq > z/ r may not be true (depending on the value of the spin parameter), as the 
ISCO may be determined by the orbital stability along the vertical direction. In this case, 
the resonances n = 1 and 2 may be exited and indeed be stronger than the resonance n = 3. 
For n = 1, the observed 3:2 ratio might be interpreted as z/u = v r + vg and z/l = v T . For 
n = 2, as v\j = 3z/ r = 3vg and z/l = 2z/ r = 2z/#. 

3.2 Forced resonances 

The equation for vertical oscillations may also include a forced resonance, in which the force 
frequency is equal to the one of the radial oscillations [46]: 

d 2 5g 

+ Qgde + [non linear terms in 5g] = h(r) cos(Cl r t) . (3-5) 
The non-linear terms allow resonant solutions for 5g, with frequencies like 

fi_ = n e - n r , (3.6) 

n + = n e + n r . (3.7) 

The observed 3:2 ratio may be explained with vg : v r = 3 : 1 (z/u = vg and z/l = vJ) or with 
vg : v r = 2 : 1 (z/u = and t'L = z^e)- I n non-Kerr backgrounds with vertically unstable 
ISCO, we have also the possibility that v r > vg and therefore, at least in principle, resonances 
like vg : v r = 1 : 2 and vg : v r = 1 : 3 may exist. 

3.3 Kepler ian resonances 

The possibility of a coupling between Keplerian and radial epicyclic frequencies might exist, 
even if it seems to be less theoretically motivated than the one in which the coupling is 
between the two epicyclic oscillations. The simplest combinations are: z^k : v r = 3 : 2 
(z^u = an d vl = v r ), UYi : f r = 3 : 1 (z^u = U K an d z^l = 2f r ), and z^k : v r = 2 : 1 (z^tj = 3i^ r 
and z^l = i^k)- For non-Kerr backgrounds, there are no new phenomena with respect to the 
standard Kerr framework, as i^k > v r is still true. 



4 The observed twin peak high-frequency QPOs in BH candidates 

In four stellar-mass BH candidates we observe two high-frequency QPOs in the X-ray flux 
and the ratio between the upper and the lower frequency turns out to be z^j : z/l = 3 : 2. For 
three of these objects, we have also a dynamical measurement of the mass, which allows us to 
use the observed twin peak high-frequency QPOs to test their nature (supposing we have the 
exact model for the production of QPOs). These three objects are listed in Tab. 1, with their 
mass M and the observed upper and lower high-frequency QPOs z;u and z/l- Regardless of 
the specific resonance model and the microphysics responsible for it, the resonance paradigm 
requires that the upper and the lower frequencies have the form 

z/U = m\v r + m 2 vg , (4.1) 
z/L = n\v r + n 2 z/0 . (4.2) 

where mi, 1712, m, and ri2 are integer (and likely as small as possible) numbers. For Keplerian 
resonance models, z/tj and z/l have clearly the form u\j = m\VYL + ' m 2Vr and z/l = niz/K + ^^r- 
Assuming the JP background with deformation parameter 63 , we can compare the theoretical 
predictions with the observed z/u and z/l for any BH candidate and define an allowed region 
on the plane spin parameter-deformation parameter for any specific set of {mi, 7712, n\, 77-2}. 
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BH binary 


M/M Q 


z/tj/Hz 


^l/Hz 


References 


GRO J1655-40 


6.30 ±0.27 


450 ±3 


300 ±5 


[55] 


XTE J1550-564 


9.1 ±0.6 


276 ±3 


184 ±5 


[56] 


GRS 1915+105 


14.0 ±4.4 


168 ±3 


113 ±5 


[36] 



Table 1. Stellar-mass BH candidates in binary systems with a measurement of the mass and two 
observed high-frequency QPOs. 

4.1 GRO J1655-40 

GRO J1655-40 is the object with the best measurements for both the mass of the BH candi- 
date and the two high-frequency QPOs, see Tab. 1. If we assume v\j and fL are the result of 
the coupling between the radial and the vertical epicyclic oscillations, we can compute v r and 
vq from Eqs. (2.13) and (2.14), consider a particular choice of {mi,m,2, ni, n-2}, and eventually 
find the allowed region on the plane spin parameter-deformation parameter on the base of 
the observed z^u and ui,. Considering only the uncertainty on the mass M and neglecting the 
one in i/tj and i-^, the allowed region for some different choices of {mi,m2,Jii,ii2} is shown 
in Figs. 3 and 4 (solid red curves) 3 . Fig. 5 shows the radial profiles of v r and vg for two JP 
space-times with €3 = 8.0, one with a* = 0.5 and the other one with a* = 0.2. In the former 
case, the condition vg : v r = 3 : 2 is satisfied at two radii and even the case ug : v T = 2 : 3 is 
possible. In the second example, we have vg : v r = 3 : 2 only at a single value of the radial 
coordinate and no vg : v r = 2 : 3. 

As we can see, there is a degeneracy between the spin parameter a* and the deformation 
parameter €3. This is the typical situation we find when we relax the Kerr BH hypothesis 
and we allow for non-vanishing deviations from the Kerr solution. In order to break this 
degeneracy, we need to use an independent constraint on a* and €3. For instance, we can use 
the one coming from the analysis of the thermal spectrum of the accretion disk of GRO J1655- 
40 during the soft-high state. This technique is called continuum-fitting method [59, 60], 
its validity is supported by a number of observational and theoretical studies [61], and its 
extension to non-Kerr background is discussed in [17, 21, 22]. In Figs. 3 and 4, the region 
allowed by the continuum- fitting method is delimited by the dashed green line 4 . 

We can repeat the same exercise to the cases with coupling between Keplerian and 
radial epicyclic frequencies. The simplest combinations are : v r = 3 : 2, 3 : 1, and 2:1. 
Fig. 6 shows the allowed regions in the plane spin parameter-deformation parameter, with 
the green dashed line still representing the bound coming from the continuum-fitting method. 

4.2 XTE J1550-564 

For GRO J1655-40, the resonance models considered in the previous subsection and the 
continuum-fitting method seem to provide consistent results only for vg : v r = 3 : 1 and z^k : 

3 When we constrain a* and 63, we cannot restrict our study to the region a*| < 1, as this bound is justified 
only for Kerr BHs [6, 7]. In the case of non-Kerr objects, the issue of the instability depends on the details 
of the structure of the object and on the gravity theory, while the creation of a very compact object with 
|a»| > 1 is possible (in the specific case of the JP space-time, at least for £3 < 0), as shown in [34, 57, 58]. 

4 The constraints coming from the continuum-fitting method for GRO J1655-40, XTE J1550-564, and 
GRS 1915+105 and shown in this paper are taken from Ref. [21]. 
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Figure 3. Constraints on the spin parameter a* and the JP deformation parameter €3 for the BH 
candidate in the binary system GRO J1655-40 from resonance models of high-frequency QPOs (red 
solid curves) and the contimium-htting method (green dashed curves). Top left panel: resonance 
vg : v r = 3 : 2 (z/y = vg and z/l = v r ). Top right panel: resonance vg : v r = 2 : 3 (i>u = v r and 
= vg). Central left panel: resonance vg : v r = 3 : 1 (z/y = Ug and 1% = vg — v r ). Central right panel: 
resonance Vg : v r = 1 : 3 [y\j = v r and = v r — vg). Bottom left panel: resonance vg : v r = 2 : 1 
(i/tj = vg + v r and = vg). Bottom right panel: resonance vg : v r = 1 : 2 (z/y = v r + vg and z^l = v r )- 
The resonances vg : v r — 2 : 3, 1 : 3, and 1 : 2 (right panels) are not allowed in a Kerr space-time. 



v r = 3 : 1 and for a non-vanishing deformation parameter £3 ~ 5 — 18. For XTE J1550-564, 
in addition to these two possibilities, even the resonances vg : v r = 2 : 1 and 1^ : v r = 2 : 1 
would be in agreement with the results of the continuum-fitting method. However, in all the 
four cases the situation is like the one shown in Fig. 7 and the degeneracy between a* and 
€3 is not broken. 
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Figure 4. As in Fig. 3, for the resonance vg : v r = 1 : 1, which is not allowed around a Kerr BH. 
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Figure 5. Profiles of the radial epicyclic frequency v r and of the vertical epicyclic frequency vg in 
the JP space-time with deformation parameter e 3 = 8.0. Left panel: spin parameter a* = 0.5; the 
resonance vg : v r — 3 : 2 is present at the radii r/M = 2.30 and 2.72, while the resonance vg : v T = 2 : 3 
is at r/M — 2.12. Right panel: spin parameter a* = 0.2; there is only one resonance vg : v r = 3 : 2, 
which takes place at the radius r/M = 5.64, and there is no resonance vg : v r = 2 : 3. 



4.3 GRS 1915+105 

For GRS 1915+105, the analysis of the thermal spectrum of the disk in the high-soft state 
in the Kerr background requires a quite high value of the spin parameter, a* > 0.97 [62], 
which makes this objects very interesting 5 . If we consider the two resonances suggested by 
GRO J1655-40 (i.e. vq : v T = 3 : 1 and v& : v T = 3 : 1), we get the two panels in Fig. 8. The 
resonance vq : v r = 3 : 1 (right panel in Fig. 8) may still produce consistent results between 
the two approaches, while the resonance z/k : v r = 3 : 1 may be ruled out. 

Fig. 9 shows instead the resonance vq : v r = 1 : 2, which is allowed only for some values 
of a* and e% (see the right panel). In this case, QPOs and continuum- fitting method provide 
consistent measurements for non- vanishing €3. That is true even for other resonances with 
v r > vq, like vq : v r = 1 : 3 and vq : v r = 1 : 1 (not shown here). In a non-Kerr background 
we may thus have several resonances that exist only if the compact object rotates sufficiently 

5 As for GRO J1655-40 and XTE J1550-564, even in the case of GRS 1915+105 I am considering the 
measurement obtained by the Harvard group. However, for GRS 1915+105 other groups have found different 
values, see e.g. Ref. [63]. 
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Figure 6. As in Fig. 3, in the case of Keplerian resonances: vk '■ v r — 3 : 1 (top panel), i^k ■ v r — 3 : 1 
(bottom left panel), and z^k : v r = 2 : 1 (bottom right panel). 




Figure 7. Constraints on the spin parameter a* and the JP deformation parameter €3 for the BH 
candidate in the binary system XTE J1550-564 from resonance models of high-frequency QPOs (red 
solid curves) and the continuum- fitting method (green dashed curves). Resonances ug : v r = 3 : 1 
(left panel) and : i> r = 3 : 1 (right panel). 



fast. 




Figure 8. Constraints on the spin parameter a* and the JP deformation parameter €3 for the BH 
candidate in the binary system GRS 1915+105 from resonance models of high-frequency QPOs (red 
solid curves) and the continuum-fitting method (green dashed curves): resonances i/$ : v r = 3 : 1 (left 
panel) and i/r '■ v r = 3 : 1 (right panel). 




Figure 9. Left panel: as in Fig. 8, for the resonance vg : v r = 1 : 2. Right panel: boundary 
delimiting the JP space-times with and without the resonance ve : v r — 1 : 2; for £3 < 0, the 
resonance vg : v r = 1 : 2 is never possible, while for e 3 > the resonance exists only for some values 
of the spin parameter (we notice that the red curve reaches a minimum at a* = 1, where 63 — > + , 
and then goes up very slowly). 



5 Discussion and conclusions 



The high-frequency QPOs observed in some stellar-mass BH candidates may be used to test 
the Kerr-nature of these objects and to confirm, or rule out, the Kerr BH hypothesis. These 
frequencies are indeed almost constant, suggesting that their value is not determined by the 
properties of the fluid accretion flow, but by the geometry of the space-time. If we assume 
the stellar-mass BH candidates are the Kerr BHs predicted by General Relativity, the high- 
frequencies QPOs may be used to estimate the spin parameter. However, at least in the 
framework of the resonance models, the estimates of the spin parameter inferred from the 
observed QPOs of three micro-quasars seem not to be consistent with the measurements 
obtained from the analysis of the thermal spectrum of the disk of the same objects, in the 
sense it is not possible to explain the observations of all the three micro-quasar by invoking a 
unique mechanism for the production of the high-frequencies QPOs (actually, this assertion 
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remains true even in all the other models proposed so far in the literature). We have thus four 
possibilities: i) the resonance models are wrong, ii) the continuum-fitting method does not 
provide reliable estimate of a*, Hi) both techniques do not work correctly, iv) these objects 
are not Kerr BHs. 

While systematic effects and/or wrong models are surely the most likely possibility to 
explain the disagreement between the two approaches, in this paper I explored the scenario 
iv) and I found two (speculative) explanations: 

1. All the observations can be explained with the resonance vq : v r = 3 : 1. As we can see 
from the central left panel in Fig. 3 and the left panels in Figs. 7 and 8, the resonance 
vq : v r = 3 : 1 may explain the high-frequency QPOs of the three micro-quasars in 
agreement with the measurements obtained from the continuum-fitting method, at the 
price of a deformation parameter €3 ~ 5 — 15; that is, these BH candidates should not 
be the Kerr BH predicted by General Relativity. GRO J1655-40 and XTE J1550-564 
would be slow-rotating compact objects and the profile of the radial and the vertical 
epicyclic frequencies would be like the one in the right panel of Fig. 5; that is, the 
resonance vq : v T = 3 : 1 would be present at a unique radius r%-\. GRS 1915+105 
would instead be a "fast-rotating" object with a* ~ 0.5 and the situation would be 
like the one in the left panel of Fig. 5: the resonance vq : v r = 3 : 1 is possible at two 
different radii, but we would observe the one occurring at the smaller one, which should 
indeed be expected to be stronger. Let us notice that the value a* « 0.5 might be close 
to the maximum value of the spin parameter for a compact object with £3 ~ 5 — 15, 
see Fig. 1 in Ref. [34]. For GRS 1915+105, the allowed region delimited by the green 
dashed line is likely in large part unphysical (equivalent to a Kerr metric with | a* | > 1 , 
i.e. it may be impossible to create such a fast-rotating objects) and therefore it is not so 
strange that the overlap between the bounds from the QPOs and the continuum-fitting 
method is a thin area near the boundary of the region allowed by the continuum-fitting 
method. 

2. The observed high-frequency QPOs are produced by different resonances determined by 
value of the the spin parameter of the compact object. Unlike in the Kerr background, for 
non-Kerr metrics such a scenario is much more appealing because there are excitation 
modes possible only if the object has a spin parameter exceeding a critical value. So, 
for a given deformation parameter, the resonance responsible for the production of the 
high-frequency QPOs would be determined by the specific value of the spin parameter 
of the compact object. This possibility can be better understood with the following 
example. In the model briefly reviewed in Subsection 3.1, resonances are possible only 
for v r jvQ = 2/n with n integer. When vq > v r , the minimum n is 3. However, for €3 > 0, 
the minimum value of n depends on the spin parameter a* . For small value of the spin 
parameter, n m [ n = 3, as in Kerr. Above a critical value of the spin parameter, which 
depends on €3, n m j n = 1 or 2. So, we may imagine a model in which GRO J1655-40 
and XTE J1550-564 present the resonance ug : v T = 3 : 1, because their BH candidates 
are slow-rotating objects and vq > v r . On the other hand, for GRS 1915+105 even 
the resonances with v r > i/g could be excited, and actually be stronger. A deformation 
parameter £3 ~ 5 — 15 is still required, i.e. the three micro-quasars should be objects 
more prolate than Kerr BHs. 

The high-frequency QPOs are potentially an excellent tool to investigate the geometry of 
the space-time around stellar-mass BH candidates and to measure the fundamental properties 



- 14 - 



of these objects. For the time being, however, we do not know the exact origin of these QPOs 
and therefore any analysis depends on the assumed model. In this paper, I have considered 
the resonance models first proposed by Abramowicz and Kluzniak, because they have a 
number of nice features, and compared the results with the measurements of the thermal 
spectrum of the disk performed by the Harvard group. All the observations of the three 
micro-quasars for which we can do this test can be explained if we admit that these objects 
are not Kerr BHs. Within the theoretical framework proposed by Johannsen and Psaltis 
with a single deformation parameter, one finds that current data demand 63 ~ 5 — 15, while 
a Kerr BHs would require £3 = 0, and for €3 > (< 0) the gravitational force on the 
equatorial plane is weaker (stronger) than the one around a Kerr BH with the same mass 
and spin. At a speculative level, we can notice that this value of 63 is consistent with another 
non-zero measurement, £3 ~ 7.5, found in Ref. [22]. The idea behind the work of Ref. [22] 
is that the steady jets observed in some micro-quasars may be powered by the Blandford- 
Znajek mechanism [64], which is one of the most appealing scenarios to explain this kind of 
phenomena. If this is the case, one should expect a correlation between the power of steady 
jets and the value of the spin parameter of the micro-quasars. Comparing current estimates of 
the jet power with the measurements of the spin parameter obtained by the continuum-fitting 
method, one finds there is no correlation if the micro-quasars are Kerr BHs, i.e. if €3 = 0, 
but a correlation is possible, and actually close to the theoretical predictions _Pj e t ~ o*, for 
e 3 « 7.5 [22]. 
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